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ABSTRACT 

Context. The Cygnus region hosts a giant molecular-cloud complex which actively forms massive stars. Interactions of cosmic rays with inter- 
stellar gas and radiation fields make it shine at y-ray energies. Several y-ray pulsars and other energetic .sources are seen in this direction. 

Aims. In this paper we analyse the y-ray emission measured by the Fermi Large Area Telescope in the energy range from 100 MeV to 100 GeV in 
order to probe the gas and cosmic-ray content over the scale of the whole Cygnus complex. The y-ray emission on the scale of the central massive 
stellar clusters and from individual sources is addressed elsewhere. 

Methods. The signal from bright pulsars is largely reduced by selecting photons in their off-pulse phase intervals. We compare the diffuse y-ray 
emission with interstellar gas maps derived from radio/mm-wave lines and visual extinction data, and a global model of the region, including other 
pulsars and y-ray sources, is sought. 

Results. The integral H i emissivity above 100 MeV averaged over the whole Cygnus complex amounts to [2.06 ±0.1 1 (stat.) (syst.)] x 10““ 
photons s“' sr“‘ H-atom"' , where the systematic error is dominated by the uncertainty on the Hi opacity to calculate its column densities. The 
integral emissivity and its spectral energy distribution are both consistent within the systematics with LAT measurements in the interstellar .space 
near the solar system. The average Xco = V(H 2 )/Wco ratio is found to be (1.68 ± 0.05 (stat.) opacity)] x 10“ molecules cm“' (K 

km s“')“*, consistent with other LAT measurements in the Local Arm. We detect significant y-ray emission from dark neutral gas for a mass 
corresponding to ~ 40% of that traced by CO. The total interstellar mass in the Cygnus complex inferred from its y-ray emission amounts to 
8 X 1 0®MO at a distance of 1 .4 kpc. 

Conclusions. Despite the conspicuous star formation activity and large masses of the interstellar clouds, the cosmic-ray population in the Cygnus 
complex averaged over a few hundred par.secs is similar to that of the local interstellar space. 

Key words. ISM: abundances - ISM: clouds - cosmic rays - Gamma rays: ISM 


1. Introduction 

Region.s w'ith conspicuous star formation activity are of great 
interest to understand the life cycle of interstellar matter and 
the properties of cosmic rays (CRs) in the Galaxy, Interstellar 
y-ray emission produced by CR interactions with the interstel- 


lar gas via nucleon-nucleon inelastic collisions and electron 
Bremmstrahlung can be used to probe their CR and gas content. 

High-energy y-ray observations have entered a new era since 
the launch of the Fenni Gamma-ray Space Telescope in 2008. 
The Fermi Large Area Telescope (LAT; Atwood et al. 2009) has 
already measured strong y-ray emission toward the 30 Doradus 




starburst region in the Large Magellanic Cloud (Abdoetal. 
2010d), and pointed out a global correlation between the y- 
ray luminosity and star-formation rate in a few normal galaxies 
(Abdo et al. 2010c). 

A primary observational target for Fermi in our Galaxy is the 
Cygnus X star-forming region, due to its proximity (~ 1 .4 kpc; 
Hanson 2003; Negueruela et al. 2008) and the availability of nu- 
merous multiwavelength observations. Named after the strong 
emission at X-ray wavelengths (Cash et al. 1980), Cygnus X is 
located around the Galactic longitude I = 80°, tangent to the 
Local Spur. It contains numerous Hit regions and OB associa- 
tions (LFyaniker et al. 2001; Le Duigou & Knodlseder 2002). It 
has long been debated whether it represents a coherent complex 
or the alignment of different structures along the line of sight. 
Recent high-resolution observations by Schneider et al. (2006) 
and Royetal. (2011) pointed out that most of the molecular 
clouds in the Cygnus X region are connected and partly show 
evidence for interactions with the massive stellar cluster Cygnus 
OB2 and other OB associations in the region . Foreground molec- 
ular clouds from the Great Cygnus Rift, at 0.6-D.8 kpc, con- 
tribute little to the large mass seen in interaction with the Cygnus 
X region itself, at 1 .4 kpc. Therefore, the molecular cloud com- 
plex appears as one of the most massive in the Galaxy. Atomic 
gas seen in these directions is probably more largely spread 
along the line of sight. 

Abdo et al. (2007, 2008) analysed Milagro measurements at 
energies > 10 TeV and reported an excess of diffuse y-ray emis- 
sion with respect to predictions based on CR spectra equivalent 
to those near the Earth; they attributed the excess to the possible 
presence of freshly-accelerated particles. 

The escape of CRs from their sources and the early propa- 
gation in the surrounding medium were so far poorly cotrained 
by observations. In particular, particles accelerated in regions of 
massive-star formation are likely to be significantly influenced 
by the turbulent environment. It is therefore interesting to inves- 
tigate how the CR populations on the scale of the massive stellar 
clusters and on the larger scale of the parent interstellar complex 
compare to each other and to the average CR population of the 
Local Spur. 

In this paper we report a global analysis of the y-ray emis- 
sion from the Cygnus region measured by the L AT in the energy 
range between 100 MeV and 100 GeV. We focus on the large- 
scale properties of the interstellar emission to probe the CR pop- 
ulation and to complement gas and dust observations at other 
wavelengths to constrain the amount of gas in different phases 
over the whole Cygnus complex. We also build an improved in- 
terstellar background framework for the study of individual y- 
ray sources which will be treated in companion papers. We dis- 
cuss interstellar emission in the star-forming region of Cygnus 
X in a dedicated paper {Fermi LAT collaboration, submitted). 


2. Data 

2.1. Gamma-ray data 

2.1 .1 , Observations and data selection 

The LAT is a pair-tracking telescope detecting photons from 
20 MeV to more than 300 GeV. The instrument is described in 
Atwood et al. (2009) and its on-orbit calibration in ,4bdo et al. 
(2009a). The LAT operates for most of the time in continuous 
sky-survey mode. We accumulated data for our region of inter- 


est from August 5 2008 (MET' 239587201) to August 5 2010 
(MET 302659202). 

We selected data according to the tightest available back- 
ground rejection criteria, corresponding to the Pass 6 Dataclean 
event class (Abdo et al, 2010e)^. In order to limit the contami- 
nation from the Earth atmospheric y-ray emission, we selected 
events with measured arrival directions within 1 00° of the local 
zenith and within 65° of the instmment boresight, taken during 
periods when the LAT rocking angle was less than 52° . 

The angular resolution of the LAT strongly depends on the 
photon energy, improving as the energy increases (Atwood et al. 
2009). Confusion at low energies is an issue since we aim 
to spatially separate the different components in the crowded 
Cygnus X region. We therefore accepted below 1 GeV only pho- 
tons which produced electron-positron pairs in the thin converter 
planes of the front section of the tracker, which provides a supe- 
rior angular resolution (Atwood et al. 2009). Above 1 GeV we 
kept all events which converted either in the front or back sec- 
tion of the tracker. 

We analysed data at Galactic longitudes 72° < / < 88° and 
latitudes -15° < b < +15°. The longitude window contains 
the interstellar complexes associated with Cygnus X; the lati- 
tude window is large enough to allow a reliable separation of 
the large-scale emission from atomic gas, isotropic background 
and inverse-Compton (IC) scattering of low-energy radiation 
fields by CR electrons. We analysed the data in the 100 MeV- 
100 GeV energy band. Below 100 MeV the instrumental system- 
atics are large (Rando et al. 2009) and the angular resolution is 
poor, whereas above 100 GeV we are limited by the low photon 
statistics, 

2.1 ,2. Removal of bright pulsars 

Three bright pulsars dominate the y-ray emission from the re- 
gion below a few GeV: the radio pulsar J2021 +3651 (Abdo et al. 
2009e) and the two LAT-discovered pulsars J202 1+4026 and 
J2032+4127 (Abdo et al. 2009b). To increase the sensitivity to 
faint sources and to the spatial structure of the diffuse emission, 
we reduced their contribution by excluding the periodic time in- 
tervals when their pulsed emission is most intense. Removing 
the intense pulsed flux helps to reduce the impact of any incor- 
rect modeling of such bright sources on the results. 

To assign pulse phases for each of the three pulsars, we pro- 
duced timing models using Tempo2 (Hobbs et al. 2006) accord- 
ing to the method described in Ray et al. (201 1)^. Fig, 1 shows 
the three light curves and the phase intervals with bright pulsed 
emission. The phase boundaries are reported in Table 1 together 
with the fraction of time in the off-pulse interval suitable for our 
study. There is a considerable level of off-pulse emission toward 
PSR J202 1+4026 that cannot be removed (Abdo et al. 2010b); 
however, given the brightness of the source, the removal of the 
on-pulse interval is useful for our aims. 

A total count-map in the off-pulse phase intervals of the three 
bright pulsars is provided for illustration in Fig. 2. To remove the 
pulsar signal without excessively sacrificing the photon statistics 
in other directions, we restricted the timing selection to a circular 
region around the pulsar position, namely to pixels in our angular 

' Fermi Mis.sion Elapsed Time, i.e, second.s since 2001 January I at 
00:00:00 UTC. 

^ Performance figures for the Dataclean event selection are given in 
the reference. 

■* For the three pulsars , the RMS of the timing residuals is below 1.1% 
of their rotational period. 






Fig. 1. Light curves and off-pulse phase intervals for the three bright pulsars. The light curves are constructed for illustration purposes 
with photons recorded in a circular region of radius 0.5° around the pulsar positions and energies > 200 MeV. 


Table 1. Off-pulse phase intervals and time fractions 


of the three bright pulsars. 


PSR 

phase interval 

time fraction (%) 

J202 1+3651 

0-0.12, 0.22-D.59, 0.7-1 

79 

12021+4026 

0.18-0.5,0.62-0.94 

64 

J2032+4127 

0-0.14,0.2-0.59,0.79-1 

74 


grid (described later in § 3 .2) the centroids of which lie within 
the energy-dependent radius 


^ cur(F) — 2 ■ 



© 0.07° 


( 1 ) 


where the symbol © indicates addition in quadrature. This is an 
approximate representation of the LAT 95% containment angle 
as a function of energy. Let us note that the accurate parametriza- 
tion of the LAT Point Spread Function (PSF) depends on energy, 
pair-conversion point in the tracker and, to a lower extent, on the 
incidence angle. The PSF is best represented by the LAT instru- 
ment response functions (IRFs), which are used later for the like- 
lihood analysis. The above acceptance-averaged approximation 
for the containment angle is only useful to calculate the radius 
fcut, and we verified that the results are insensitive to reasonable 
variations of this parameter. 

To take into account the cut on pulsar phases, for each di- 
rection in the sky and energy the exposure (see again § 3.2) was 
multiplied by the remaining livetime fraction. The remainder of 
the pulsar emission was included in the model using: 


a point source to represent emission in the off-pulse interval; 


- a second point source, for which the number of expected 
counts is set to null at r < rcut{E) from the pulsar position, to 
represent on-pulse y-rays spilling over at r > r,, „,(£■) in the 
tails of the PSF. 

The two sources have free independent fluxes in each energy 
bin of the analysis. This is particularly important to account for 
the different spectra of the on-pulse and off-pulse y-ray emission 
and also to compensate for any mismatch between the tails of the 
model PSF and the emission at large angles from the brightest 
sources in the region. 

Since the three pulsars have exponential spectral cutoffs near 
2-3 GeV (Abdo et al. 2010b) the phase selection was not ap- 
plied above 10 GeV where the level of pulsed emission is low 
and each pulsar was accounted for by a single point source. On 
the other hand, given the large statistics but poor angular reso- 
lution at low energies (more than half of the region of interest 
would be subject to on-pulse event removal) we selected off- 
pulse photons for the whole region below 316 MeV In this 
case no “spill-over” source was necessary. 

2.2. Ancillary data 

2.2.1. Radio/mm-wave lines: neutral gas 

Neutral atomic hydrogen, H i, is traced thanks to its 21 -cm line. 
Where available* we used data from the Canadian Galactic Plane 
Survey (CGPS; Taylor et al. 2003) rebinned onto the 0.125° x 
0.125° grid used for the other maps. Elsewhere, we used data 

See § 3.2 for the definition of the energy grid used in the analy.sis. 

° The CGPS coverage is almost complete at -3.5° < b < +5.5° for 
this longitude range. 
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Fig. 2. Total count map in the energy range 100 MeV-100 GeV, 
binned over a 0.125° x 0.125° grid in Galactic coordinates in 
Cartesian projection. Data were selected according to the crite- 
ria described in the text (§ 2.1.1) and in the off-pulse phase in- 
tervals of the three bright pulsars (§ 2.1 .2), whose positions are 
marked by diamonds. Counts are saturated between 0 and 50, 
and smoothed for display with a Gaussian kernel of cr = 0.25°. 


from the Leiden/ Argentine/Bonn (LAB; Kalberla et al. 2005) 
survey, with a coarser binning of 0,5°. We checked the consis- 
tency of the two survey calibrations in the overlap region. 

Molecular hydrogen cannot be observed directly in its most 
abundant cold phase. The velocity-integrated brightness temper- 
ature of the '-CO 2.6-mm line, Wco, is often assumed to lin- 
early scale with the /V(H 2 ) column density. We used CO data 
from the composite survey by Dame et al. (2001). filtered with 
the moment-masking technique (Dame 201 1) in order to reduce 
the noise while preserving the faint cloud edges and keeping the 
resolution of the original data. 

The Doppler shift of radio/mm-wave lines can be used to 
kinematically separate the Cygnus complex from two faint seg- 
ments of the Perseus and outer spiral arms seen beyond Cygnus 
in the same direction. We applied the kinematic separation pro- 


cedure illustrated by Abdo et al. (201 Of), starting from a prelim- 
inary boundary located at a Galactocentric radius* R = 9.4 kpc 
and then adapting the separation to the cloud structures and 
correcting for the spill-over due to the broad velocity disper- 
sion of H I lines. The separation into two regions is sufficient to 
model the interstellar y-ray emission in Cygnus since Abdo et al, 
(20100 and Ackeimann et al. (201 1) did not find significant gra- 
dients of the gas y-ray emissivities in the outer region of the 
Milky Way, We applied the kinematic separation procedure to 
prepare maps of the column densities of atomic gas, A(H i), and 
of Wco ■ The maps are shown in Fig . 3 for H i and Fig . 4 for CO , 
Note that all the gas maps mentioned along the paper have > 10° 
borders around the analysis region used to properly convolve the 
model with the LAT PSF. 

Substantial uncertainties in the determination of N (H i) arise 
from the choice of spin temperature for the optical depth correc- 
tion. We adopted as baseline case a uniform Ts = 250 K, which 
is the average spin temperature that best reproduces the blend- 
ing of cold and warm atomic gas according to observations of 
emission-absorption H I pairs in the region covered by the CGPS 
(Dickey et al. 2009). Other values 100 K<Ts <oo will be con- 
sidered later to evaluate the related systematic uncertainties af- 
fecting the results of our analysis. 

2.2.2. Visual extinction: dark neutral gas 

Multiwavelength observations indicate that the combination of 
the Hi and CO lines does not properly trace the total col- 
umn densities of the neutral interstellar medium (ISM) (e.g. 
Magnanietal. 2003; Grenier etal. 2005; Abdoetal. 2010f; 
Langeretal. 2010; Adeetal, 2011). Since the work by 
Grenier et al. (2005), dust tracers have been used in y-ray anal- 
yses to complement the Hi and CO lines, under the assumption 
that dust grains are well mixed with gas in the warm and cold 
phases of the ISM and therefore provide an estimate of total gas 
column densities. Grenier et al, (2005) and Abdo et al. (20101) 
adopted the E(B - V) color excess map by Schlegel et al. (1998) 
as a tracer of the total column densities, and used the E(B - V) 
residuals -i.e. E(B - V) minus the best-fit linear combination of 
N(H l) and Wco maps- as a tracer of the dark-gas column densi- 
ties in nearby clouds of the Gould Belt. 

The use of the E(B - V) map is problematic in the Cygnus X 
region for two reasons: 

- numerous infrared point sources contaminate the map; 

- the temperature correction used by Schlegel et al, (1998) to 
derive the dust column-density map from IRASflSSA mea- 
surements is highly uncertain in regions of massive-star for- 
mation because of the enhanced radiation fields. 

We have therefore adopted the visual extinction Ay as 
derived from the reddening of near-infrared sources in the 
2MASS catalog (Skrutskie et al. 2006). The Av maps produced 
by Rowles & Froebrich (2009) and Froebrich & Rowles (20 1 0) 
were used for Ay < 5 mag. They exhibit saturation effects 
at higher extinction values, so we complemented them with 
an Av map obtained from 2MASS data using the code and 
method developed by Schneider et al. (2011). The latter use 
the Besanfon stellar population model (Robin & Creze 1986; 
Robin et al, 2003) to filter out the contribution from the fore- 


* Based on the assumption of a flat rotation curve with solar radius 
RQ = 8.5 kpc and Galactic rotation velocity at the solar circle Vo = 
220 km s-‘ 





Fig. 3. Maps of iV(Hi) column densities in the Cygnus complex in the Local Spur (left) and in the outer Galaxy (right), under the 
assumption of a uniform spin temperature of 250 K. The colour scales with A?(Hi) in units of 10^° atoms cm~^. The maps were 
smoothed with a Gaussian kernel of cr = 0.25° for display. 


ground bluest stars’^. The second Ay map was built in a 12° 
region centered on (/, b) = (80°, 0°); compared with the first 
set of maps, it presented an offset of ~ 0,46 mag at low 
extinction. We constructed the final Ay map from the direct 
Rowles & Froebrich (2009) data below 5 mag and from the sec- 
ond map, offset by 0.46 mag, at higher extinction. 

The Ay map was binned onto the same 0. 1 25° x 0. 1 25° grid 
in Cartesian projection as the other maps. The Ay map was fit- 
ted with a linear combination of the N(H l) and Wco maps pre- 
viously described. The input Ay map minus the best-fit linear 
combination of the A(H l) and Wco maps yielded the Ay excess 
map, Ay_exc. which will be used to trace the dark neutral gas. 
Only residuals coiresponding to input Ay > 0.3 mag were kept 
to limit the noise off the plane. The Ay excess map is shown in 
Fig. 5. 


^ To do so, a distance from the observer needs to be assumed for the 
clouds under consideration. We verified that variations of the order of 
a few hundred parsecs do not significantly change the results presented 
in the paper. 


2.2.3. Microwave emission: ionized gas 

Away from H II regions around massive stars and stellar clus- 
ters, the ionized gas constitutes a layer of characteristic height 
> 1 kpc over the Galactic plane with little mass compared to 
the neutral phases (Cordes & Lazio 2002). Therefore, it has of- 
ten been neglected in previous y-ray studies. However, we find 
in the Cygnus X region many conspicuous H II regions excited 
by the intense radiation fields of the numerous massive stars 
(Uyaniker et al. 2001 ; Paladin! et al. 2003). 

Ionized gas masses can be traced by free-free emission fol- 
lowing the prescription by Sodroski et al. (1989, 1997) to derive 
the A'(H ll) column densities: 



where Iff is the free-free emission intensity at the frequency v, T\. 
is the electron temperature and «eff the effective electron num- 
ber density. We adopted a free-free emission map derived from 
the 7-year WMAP data in the Q band (40 GHz) by Gold et al. 
(201 1) using the maximum entropy method from the prior tem- 
plate given by the extinction-corrected Ha map by Finkbeiner 
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Fig. 4. Maps of VFco intensities in the Cygnus complex in the Local Spur (left) and in the outer Galaxy (right). The colour scales 
with IVco in units of K km s“* above 1 .5 K km s“' . The maps were smoothed with a Gaussian kernel of cr = 0.25° for display. 


(2003). It was rebinned onto the 0.125° x 0.125° grid used for 
the other maps, as shown in Fig. 6. 

3. Analysis 

3.1. Analysis model 

3.1.1. Diffuse emission 

Since the bulk of Galactic CRs in the relevant energy ranges are 
expected to be smoothly distributed over scales exceeding the 
typical dimensions of interstellar clouds and to penetrate uni- 
formly all the phases of the ISM, the y-ray emission produced 
by CR-gas interactions can be modelled to first order as a linear 
combination of the gas column densities summed for the differ- 
ent phases and different regions along the line of sight. 

Ionized gas, with a total mass over the region of 0.4 x 
10® (ncff/1 cm-^)-‘ Mo for Te = 10^ K, represents less than 4% 
of the total atomic mass present in the Cygnus complex (assum- 
ing tieff = 2 - 10 cm~^, Sodroski et al. 1997). The corresponding 
column densides are highest in the massive-star forming region 
of Cygnus X where we detected a bright and hard e.xtended y- 
ray source powered by freshly-accelerated particles, that will be 
called hereinafter “the cocoon” {Fermi LAT collaboration, sub- 


mitted). The free-free emission map was significantly detected 
in addition to the other interstellar components, but only at the 
expense of an unusually hard spectrum. In order to model the 
entire region, we introduced an extended source to account for 
the cocoon, as described in 3.1 .2. The latter was found to pro- 
vide the best fit to the LAT data, yielding a larger maximum- 
likelihood value than the free-free emission template. Since the 
cocoon source overlaps most of the ionized clouds it absorbs 
their contribution to the y-ray emission; ionized gas was there- 
fore not included in the baseline model through the free-free 
emission template. The Cygnus X region is treated in detail in 
a companion paper {Fermi LAT collaboration, submitted). The 
results presented in this paper were checked against the inclu- 
sion of the free-free emission template in the model. 

The interstellar IC emission is produced by interactions of 
CR electrons and positrons with the low-energy interstellar ra- 
diation field (ISRF). To account for large-scale IC emission 
from the Milky Way we adopted a template calculated using the 
GALPROP CR propagation code* (Strong & Moskalenko 1998; 
Strong etal. 2007), run 54_87Xexph7S. The IC emission was 
calculated on the basis of a CR electron spectrum consistent with 


® http://ga!prop..stanford.edu 






Fig. 5. Av excess map (magnitudes) obtained from the optical 
extinction Ay estimated from 2MASS data minus the best-fit 
linear combination of the iV(H l) and IFco maps shown in Fig. 3 
and 4, respectively. The map was smoothed using a Gaussian 
kernel of cr = 0.25° for display. 


recent measurements at the Earth (Abdo et al. 2009c) and the 
new calculation of the Galactic ISRF by Porter et al. (2008). 

Local radiation fields could leave unmodelled structures in 
IC emission, notably in the massive-star forming region of 
Cygnus X (e.g. Orlando & Strong 2007), In the companion pa- 
per we show that an upper bound to the IC emission from the 
stellar and interstellar low-energy radiation fields upscattered by 
CR electrons with the local spectrum is 2 orders of magnitude 
fainter than the cocoon emission, which in turns is fainter than 
the emission from the neutral gas (Fig. 9). CR electron sources 
within Cygnus X could further enhance the IC y-ray yield. Any 
enhanced IC contribution from the inner region is accounted for 
in this analysis by the extended cocoon source and it should not 
bias the determination of the gas emissivities we aim at studying 
here. 

The diffuse emission model is completed by the isotropic 
background which combines the residual backgrounds from 


Fig. 6. Free-free emission intensities from WMAP data. The 
color scales with brightness temperature in mK. The map was 
smoothed using a Gaussian kernel of m = 0.25° for display. 

misclassified CR interactions in the LAT and the isotropic, 
presumably extragalactic, y-ray emission (studied in detail in 
Abdoetal.2010e). 

3.1.2. Sources 

We included in the model the identified sources in the region 
of interest: Cygnus X-3 (Abdo et al. 2009f), PSR J1957+5033 
(Saz Parkinson et al, 2010) and PSR J2030-I-364I (Camiloet al., 
in preparation), in addition to the three bright pulsars as dis- 
cussed in § 2.1.2. 

We also iteratively included significant 1 FGL point sources 
(Abdoetal. 2010a) either associated with Active Galactic 
Nuclei (AGN) or characterized by variability or both: the 
sources were added with decreasing brightness: J21 16. 1-1-3338, 
J200I.l-f4351, J2027 .6+3335, J21 15.5+2937, J2015.7+3708, 
J2029 .2+4924, J2012.2+4629 and J2128.0+3623. The iterative 
procedure is useful to stabilize the likelihood fitting procedure 
and to assess the significance of sources added at each step. 








We detected extended y-ray emission above the global 
interstellar emission model discussed here associated with 
the supernova remnants known as the Cygnus Loop (G74.0- 
8.5, e.g. Sunetal. 2006) and y Cygni (G78.2+2.1, e.g. 
Ladouceur & Pineault 2008) and with the inner 100 pc of the 
Cygnus X complex. They are discussed in detail elsewhere. We 
briefly summarize here how these extended sources are mod- 
elled. For each of them we have tested different models and 
therefore verified that their presence does not bias the results 
concerning the properties of large-scale interstellar emission 
presented in the paper. 

The Cygnus Loop was modelled using a ring centered at 
{l,b) - (74.1°, -8.5°) and with inner/outer radii of 0.7° and 
1.6°, respectively, which best reproduces y-ray emission from 
the Cygnus Loop (Katagiri et al. 201 1). 

We included two sources in the region of y Cygni in addition 
toPSR J202 1+4026: 

- a uniform disc centered at [I, b) = (78.2°. +2. 1°) and a radius 
of 0.5° (G78.2+2.1 ; Green 2009); 

- a 2D Gaussian corresponding to the moderately extended 
TeV source^ VER 2019+407 (Weinstein et al. 2009). 

We detected extended y-ray emission toward the inner ~ 
100 pc of Cygnus X, which is effectively treated here as a 
source named “the cocoon”. We discuss the nature of the co- 
coon and its relation with CR acceleration in the massive- 
star forming region in a dedicated paper {Fermi LAT collabo- 
ration, submitted), where we determine a Gaussian centered at 
(/, b) = (79.6° ± 0.3°, 1.4° ± 0.4°) with a cr = 2.0° ± 0.2° width 
to be the model providing the best fit to the LAT data. As noted 
above, the cocoon source effectively accounts for the contribu- 
tion from ionized gas and from enhanced IC emission in the 
Cygnus X massive-star forming region, as well as for the pres- 
ence of locally-accelerated CRs. 

The spatial distribution of the sources included in the analy- 
sis model is summarized in Fig. 7. 

3.1 .3. Summary of the analysis model 

To summarize, the y-ray intensities / (photons cm~^ s“* sr“*) are 
modelled in each energy bin by 

2 

/(/, fc) = £ [qm., ■ N(U !)(/, b), + qco., ■ Wcod, b),] + 

+ ^Av ■ Av ,exc (/, b) + IC(/, b) "f /iso + 

+ YjSjihb). (3) 

J 

The sum over i represents the combination of the two regions: 
I) Cygnus complex and 2) outer Galaxy. The free parameters 
of the diffuse emission model are the emissivities per hydrogen 
atom, ^Hi.i (s“' sr"'), the emissivities per unit of Wco intensity, 
qco, I (cm'^ s*' sr“‘ (K km s~')“‘), the emissivity per Ay ex- 
cess unit ?Av (cm'- s~' sr'‘ mag~' ) and the isotropic intensity 
/iso (cm“^ s'* sr'’). IC stands for the aforedescribed IC emis- 
sion model. The sum over j represents the combination of the 
sources, either point-like or extended, as described in § 3.1.2, 
including a free parameter (flux normalization) independently 
for each of them, 

® The 2D Gaussian fitted to TeV data by Weinstein et al. (2009) pro- 
vides a better fit to LAT data with respect to the coincident point source 
IFGL J2020+4049, 



Fig. 7. Sources included in the analysis model. Diamonds mark 
the positions of the three bright pulsars for which photons 
were phase-selected (§ 2.1.2). X points mark the positions of 
other identified sources. Crosses correspond to IFGL sources 
either associated to AGN, or variable, or both, confirmed by our 
analysis. The blue circles correspond to the rims of the tem- 
plates adopted to model SNRs, the Cygnus Loop (G74.0-8.5) 
and y Cygni (G78.2+2.1). The red circle marks the centroid 
of VER 2019+407 (whose extension is not appreciable in this 
large-scale view). The magenta circle represents the l<r contour 
of the Gaussian used to model the cocoon. 


3.2. Analysis method 

The model was fit to LAT data by using a binned maximum- 
likelihood with Poisson statistics'® independently over several 
energy bins. We used a 0.125° x 0.125° binning in Cartesian 
projection, comparable with the LAT angular resolution at the 
highest energies. We considered three energy bands: low (100 
MeV- 1 Ge V) , raid ( 1 GeV- 1 0 Ge V) and high energies ( 1 0 GeV- 
100 GeV). The low and mid-energy bands were further divided 
in four logarithmic-spaced energy bins; the higher-energy band 
in two because of the limited statistics" . 

A.S implemented in the standard LAT analysis tools 89-18-85, 

' ' The bounds of the energy bins are reported in Table 2. 




The analysis was based on the post-launch IRFs of the P6_V3 
series, which take into account efficiency losses due to pile-up 
and accidental coincidence effects in the detector (Rando et al. 
2009). 

To perform the convolution with the LAT PSF, a power-law 
spectrum with index 2.1 was assumed for the gas maps and 
other sources modelled by geometrical templates; the results do 
not significantly depend on this value. For all other sources we 
used power-law spectra with the spectral index reported in the 
IFGL Catalog. For the pulsars included in the LAT pulsar cata- 
log (Abdo et al . 20 1 Ob) we used the spectral functions described 
therein. 

4. Results and discussion 

4.1. Summary of the results and uncertainties 

The y-ray residuals corresponding to the best-fit model are 
shown in Fig. 8. They indicate that the model satisfactorily re- 
produces the morphology of the y-ray emission at angular scales 
larger than the LAT PSF in all the energy bands. Localized posi- 
tive residuals are still present. Some of them coincide with unas- 
sociated IFGL sources; others are associated with sources of the 
2FGL catalog'^, notably 2FGL J2018.0 4- 3626 also coincident 
with the TeV source MGRO J20 19+27 (Abdo et al. 2007). We 
verified that including sources accounting for those residuals in 
the analysis model would not significantly affect the determina- 
tion of y-ray emissivities associated with the different gas com- 
ponents summarized in Table 2. 

Fig. 9 shows the y-ray spectral energy distribution measured 
by the LAT over the whole region of interest. LAT measurements 
are compared with the final model, and the different components 
are outlined. The data sample is dominated by emission from in- 
terstellar gas in the Cygnus complex. The largest contributor is 
Hi. Emission associated with CO and Ay excesses exceeds the 
signals from individual sources for the whole energy range con- 
sidered. The cocoon has a very hard spectrum and becomes com- 
parable to emission from CO-bright gas at energies > 10 GeV. 

All the results presented so far are based on the assumption 
of a uniform Hi spin temperature of 250 K. To gauge the im- 
pact of the optical depth correction of Hi data on the results, 
we repeated the analysis with other assumptions, Ts = 400 K is 
considered since it is the value best reproducing pairs of emis- 
sion/absorption H I spectra over most of the regions analysed by 
Dickey etal. (2009), although they found that Ts = 250 K is 
preferred in the region covered by CGPS data. Ts = 125 K 
is considered as it has long been used for y-ray studies (e.g. 
Bloemen et al. 1984; Hunter et al. 1997; Strong et al. 2004), We 
also considered two extreme assumptions: a low'^ Ts = 100 K. 
and the optically thin approximation (equivalent to infinitely 
high spin temperature). Fig. 10 shows the maximum likelihood 
profile obtained for the final model as a function of Ts . The 
results support the average spin temperatures of a few hun- 
dred K deduced from radio absorption/emission measurements 
by Dickey etal, (2009), implying a mix of < 25% cold and 
> 75% warm Hi. 

In order to test the robustness of the results against the pres- 
ence of ionized gas beyond the extended cocoon source, we re- 
placed the latter in the baseline model with the free-free emis- 

The preparation of the 2FGL catalog ran in parallel with the anal- 
ysis reported in this paper. The source list is now available from link 
http://fermi.gsfc.na.sa,gov/ssc/data,taccess/lat/2yr„catalog/. 

The spin temperature is higher than the brightness temperature, 
measured > 100 K along many directions in the region. 
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Fig. 10. Maximum likelihood obtained as a function of the uni- 
form spin temperature adopted for the optical depth correction 
of H I data. Values are offset so that log-likelihood is null for a 
spin temperature of 100 K, 

sion template. All the results on the gas emissivities were found 
to be consistent with the values listed in Table 2 within statistical 
errors . 

The large-scale 1C model introduced in § 3.1 .1 is affected by 
considerable uncertainties related to the distribution of CR den- 
sities and of the ISRF in the Galaxy. The intensity of IC emission 
expected over our region of interest is comparable to that from 
interstellar gas. The latter, however, is highly structured and has 
a smaller characteristic height above the Galactic plane, and can 
be therefore reliably determined in the likelihood fit. We verified 
that completely neglecting the large-scale IC emission leads to 
negligible variations in the emissivities of CO-bright gas and Ay 
excesses, and to variations lower or comparable to statistical un- 
certainties for the emissivities of the atomic hydrogen, which is 
less structured than the other ISM components and has a larger 
characteristic height. 

Other systematic uncertainties are due to the LAT instrument 
response. The uncertainties in the y-ray selection efficiency are 
estimated to be 10% at 100 MeV, 5% at 560 MeV, and 20% 
above 10 GeV for the IRFs we used here (Abdo etal. 20I0e). 
The Monte Carlo-based PSF used for this study is known to not 
accurately reproduce in-flight data over the whole energy range 
considered. We verified by means of dedicated simulations that 
this does not significantly affect the determination of the gas 
emissivities considered for the discussion. The energy disper- 
sion is routinely neglected in the likelihood fitting of LAT data 
for limitations in computing power: Monte Carlo simulations in- 
dicate that this approximation causes a bias of the order of 10% 
at 100 MeV decreasing to < 5% above 200 MeV. 

4.2. H I emissivity and CR densities 

The Hi emissivity per hydrogen atom relates to the average 
CR density in each of the regions considered. LAT measure- 
ments (Abdo et al. 2009d) show that the H l emissivity spec- 
trum in the local ISM is consistent with production via electron 
Bremsstrahlung and nucleon-nucleon interactions by CRs with 
a spectrum consistent with that directly measured in the Earth 
neighborhood. 

The integrated y-ray emissivity > 100 MeV we measure in 
the Cygnus complex amounts to [2.06±0, 1 1 (stat.) (syst.)] x 
10“25 j-i gj.-i pjg ]] sjjQws the Hi emissivity spectrum ob- 
tained for the Cygnus complex and compares it with the expec- 
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Fig. 8. Residuals (data-model). Left: low energies (100 MeV-1 GeV); center: mid-energies (1 GeV-10 GeV); right: high energies 
(1 0 GeV- 100 GeV). Units are approximate standard deviations (square root of model counts) saturated between ±3cr and smoothed 
for display with a Gaussian kernel of cr = 0.25°. In each panel the blue circle in the top left comer represents the effective LAX 
PSF 68% containment circle for the event selection used in the analysis (averaged over the corresponding energy range assuming 
a power-law spectram with index 2.1). Diamonds mark the positions of bright pulsars for which phase-selection was applied, as in 
Fig. 2. Crosses mark the positions of unassociated IFGL sources coincident with positive residuals; the X point mark the position 
of 2FGLJ201 8.0 + 3626, coincident with a hot spot in the residual map 1 GeV-10 GeV. 

Table 2. Best-fit parameters describing emission from interstellar gas (Eq. 3) under the assumption of a uniform H i spin temperature 
Ts = 250 K. 


energy bin“ 

/7 ^ 

<?CO, 1 

^Hl,2 

9co,2 ° 

9av " 

0.1-0,178 

7.9+1.1 

3. 3+ 1.0 

8.5+1. 4 

0.00+0,06 

10+30 

0.178-0.316 

5.9+0.3 

2.2+0.3 

4.7+0.4 

0.000+0.002 

19+5 

0.316-0.562 

3.27+0.11 

LI 6+0 .08 

3.23+0.14 

0.00+0.06 

11.0+1.8 

0.562-1 

1.95+0.06 

0.59+0.04 

1.72+0.10 

0.5+0 .3 

6.2+0.7 

1-1,78 

0.98+0.03 

0.328+0,016 

0,74+0.04 

0,12+0.13 

2,6±0.3 

1.78-3.16 

0.389+0.016 

0.141+0.008 

0.36+0.02 

0,02+0.06 

0.86+0,15 

3.16-5.62 

0.151+0.005 

0,044+0.004 

0.113+0.013 

0.02+0.03 

0.39+0.08 

5.62-10 

0.050+0.003 

0.016+0.002 

0.046+0.006 

0.000+0,005 

0.15+0.04 

10-31,6 

0.0085+0.0015 

0.0059+0.0010 

0.021+0.003 

0.002+0.008 

0.043+0.019 

31.6-100 

0.0024+0.0007 

0.0016+0,0004 

0.0007+0.0014 

0.002+0.003 

0.002+0.007 


Notes. Subscripts refer to the two regions separated in analysis: 1 ) the Cygnus complex in the Local Spur, 2) the outer Galaxy. Some parameters 
are poorly determined but they are reported for completeness. 

M GeV 10 - 2 ’ 5-1 5J-1 M (0-6 ^-2 5-1 5^-1 5-i)-i M iq-6 ^-2 j-i 5^-! mag-' 


tations for the local interstellar spectrum estimated in Abdo et al. 
2009d. The latter is compatible (within 10%) with LAX observa- 
tions at mid latitudes in the third Galactic quadrant in the energy 
range 100 MeV-10 GeV (Abdo et al. 2009d). The spectrum is 
presented for a uniform spin temperature Tg = 250 K; system- 
atic uncertainties due to the H i opacity correction and to the y- 
ray selection efficiency are added in quadrature for display. The 
latter give a non-negligible contribution only at energies larger 
than a few GeV. 

The eraissivity of atomic gas in the Cygnus region, averaged 
over ~ 400 pc, is consistent with the local emissivity in the 


100 MeV-100 GeV energy range, except for the deviant point 
at 10 - 30 GeV. The latter can be explained by the difficulty in 
distinguishing the H i and CO components at high energies, for 
which the y-ray statistics are limited and the hard cocoon source 
is brighter (Fig. 9). The emissivity spectrum implies that the CR 
spectra in the relevant energy ranges (~ 1 - 100 GeV/n for nu- 
cleons) are similar to those measured in the vicinity of the Earth 
and inferred from y-ray observations in the nearby interstellar 
space within 1 kpc. 

The variations in average CR densities along the Local Spur 
between the dense Cygnus complex, two segments in the sec- 
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Fig. 9. Spectral energy distribution of y-ray emission measured by the LAT compared with our best-fit model. Statistical errors 
only. We separately show the different components of the interstellar emission model, point sources and extended objects. The 
curve corresponding to y Cygni combines the contributions of the off-pulse source, the disc associated with the remnant and the 2D 
Gaussian accounting for VER 20 1 9+407 . 
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Fig. 11. H I emissivity spectrum in the Cygnus complex. Points: 
the best-fit estimate for the spin temperature Ts = 250 K. 
Hatched rectangles: systematic uncertainties taking into account 
H I opacity and y-ray selection efficiency. Line: model of the lo- 
cal interstellar spectrum by Abdoetal. 2009d (with a nuclear 
enhancement factor of 1 .84, Mori 2009). 


ond and third quadrants that exhibit ~ 2 lower surface densities 
of gas (Abdo et al. 201 Of: Ackermann et al. 201 1), and the mid- 
latitude diffuse medium with a factor ~ 5 lower surface density 
(Abdoetal. 2009d) are constrained to be within 10% - 35%: 
this is difficult to reconcile with the idea of a dynamical cou- 
pling between gas and CR densities (e.g. Bertsch et al. 1993; 
Hunter et al. 1997). They are consistent on the other hand with 


the small arm-interarm emissivity contrast estimated from LAT 
data in the third Galactic quadrant (Ackermann et al. 201 1). 

In spite of the large column densities of gas, exceeding 10^^ 
atoms cm“^ over many directions within the Cygnus complex, 
we find no hints of exclusion of CRs from the densest parts of 
the atomic clouds. 

Due to the bright foreground of the Cygnus complex and 
individual sources, studying in detail the gas emissivity in 
the outer disc of the Milky Way is beyond the scope of this 
study. However, the ratio of the integrated Hi emissivity of 
the outer region over that in the Local Spur is (90 ± 7)%, in 
very good agreement with the results by Abdoetal. (201 Of) 
and Ackermann et al . (20 1 1 ) . It confirms in another direction the 
presence of large CR densities beyond the solar circle. 

Located at a distance of ~ 1.4 kpc and I = 80° , the Cygnus 
complex lies at R = 8.4 kpc from the Galactic center, at a slightly 
smaller radius than the solar system. We measure in this direc- 
tion an emissivity which is consistent with other values found 
in the Local Spur and in three outer segments of the Milky Way 
with Galactocentric radii up to ~ 15 kpc. It yields a decrease 
in H I emissivity < 60% over ~ 6 kpc in Galactocentric radius 
across and beyond the solar circle. Yet, those measurements span 
a small range in azimuth around the Galactic center and may not 
be representative for comparison with axisymmetric propagation 
models in order to study the CR gradient across the solar circle . 

4.3. CO-bright molecular gas 

If molecular and atomic gas are illuminated by the same CR 
fluxes vve expect the emissivity per hydrogen molecule to be 
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Fig. 12. Top: emissivity per VPco intensity unit versus emissivity 
per hydrogen atom in the Cygnus complex (for Ts = 250 K). 
The points correspond to the different energy bins; the emissivi- 
ties decrease with increasing energy. The red line gives the best 
linear fit taking into account uncertainties on both axes. Bottom: 
residuals in units of standard deviations as a function of H l emis- 
sivity. 


twice the emissivity per hydrogen atom and we can therefore 
calibrate the Xqo ratio. We performed a linear fit, qco.i = 
q + 2Xqo ■ qnt,i, taking into account the uncertainties on both 
emissivities to derive the best linear relation shown in Fig, 12; 
we also give the residuals in units of standard deviations. A good 
linearity is found in the 0.1-10 GeV energy range. The highest- 
energy (lowest-emissivity) points show < 3cr excess of emission 
associated with CO with respect to the Xco ratio determined at 
low energies. The high CO emissivity recorded at 10-30 GeV 
(second point) corresponds to a low emissivity in Hi (Fig. 1 1) 
and may result from a fluctuation in the difficult spatial separa- 
tion between the atomic and molecular components when pho- 
tons are sparse and in the presence of the hard 2° source which 
partially overlaps the CO peaks. Up to 10 GeV the linearity is 
good and therefore there is no sign of CR exclusion from the 
dense cores of this giant molecular complex. 

The slope of the best-fit linear relation provides a value 
of Xco = (1-68 ± 0.05) X 102» cm'^ (K km s-')-' in the 
case of Ts = 250 K. We obtain Xco = (1.58 ± 0.04) x 10^® 
cm“^ (K km s“')“' in the limit of small Hi optical depth and 
Xco = (2.55 ± 0.08) X 10^® cm“^ (K km s~')“‘ in the case of 
Ty = 100 K. The uncertainties in ^hi associated with the H i spin 
temperature are particularly severe for the high-density clouds 
of the Cygnus complex. High optical depths (low spin temper- 
atures) imply a large increase in ATHi), therefore substantially 
lower CR densities. Given the y-ray luminosity of the molecular 
clouds, this subsequently implies a significant increase in their 


estimated masses*'*. The systematic errors on the y-ray selection 
efficiency cancel out to first order in the estimate of the Xco ra- 
tio. 

The conversion factor Xco = [1-68 ± 

0.05 (stat.) lojolHl opacity)] x 10“ cm'^ (K km 8“*)“' 
is consistent with other LAT measurements in the Local Spur 
which range from 1 .5 to 2 xl0“ cm^‘ (K km s^*)^‘ (Abdo et al, 
201 Of; Ackermann et al. 2011). From the different y-ray mea- 
surements in the Galactic plane the Xco ratio at the solar circle 
appears well defined. It is, however, significantly larger than 
in nearby well-resolved clouds off the plane in Cassiopeia and 
Cepheus (Abdoetal. 2010f). Whether the discrepancy is due 
to the sampling resolution or to an intrinsic Xqo variation at 
different scales inside a cloud will be investigated in the future. 

Using the Xco ratio, we estimated the CO-bright molecular 
mass in the complex. For this purpose we considered the region 
at 74° < I < 86°, -5° < 6 < 8°, where most of the gas as- 
sociated with the Cygnus complex is located. Assuming a dis- 
tance of 1 .4 kpc and a mean atomic weight per hydrogen atom 
in the ISM of 1.36, we obtain a mass 2.3 j x 10® MO (where 
the uncertainties are dominated by the Hi opacity correction). 
This value (taking into account the different assumption on the 
distance) is consistent with the results by Schneider et al. (2006) 
based on higher-resolution, multi-isotopologue CO observations 
and it depicts Cygnus as a super-massive molecular complex. 

Due to the small amount of CO-bright molecular gas in the 
outer region of the Milky Way in this longitude window (Fig. 4), 
the determination of its emissivities is extremely sensitive to the 
details of the model (including point sources) and we do not con- 
sider it for scientific interpretation. 

4.4. Dark neutral gas 

The inclusion of the Ay excess map in the model corresponds 
to an increase of 250.6 in the logarithm of the likelihood (for 
ten additional degrees of freedom). This corresponds to a signif- 
icant detection of y-ray emission associated with Ay excesses, 
formally equivalent to a ~ 2 Im confidence level. 

Fig. 13 shows the emissivity per Ay,exc unit, ^Av> versus the 
emissivity per hydrogen atom, qm, in the Cygnus complex. A 
good linear correlation is found between the two emissivities 
over three decades in energy, proving that y-ray emission asso- 
ciated with Ay excesses originates from the same physical pro- 
cesses as that associated with H l. Ay residuals therefore trace 
interstellar gas. 

With a procedure analogous to that adopted to estimate Xco, 
we can use the emissivity per hydrogen atom to calibrate the 
dust-to-gas ratio in the dark neutral phase Xav = V(H)/Ay,exc- 
We obtain Xav = (28 ± 2) x 10“ cm^^ mag'"' in the case 
of Ts = 250 K, Xav = (48 ± 3) x 10“ cm“^ mag"' in 
the case of = 100 K and Xav = (27 ± 2) x 10“ cm"^ 

mag"' in the case of optically thin medium. Xav is therefore 
[28 ± 2 (stat.) (H i opacity)! x 10“ cm"-^ mag"' . 

Assuming a standard total-to-selective extinction ratio Ry = 
Av/E(B - V) = 3.10 (Wegner 2003), the dust-to-gas ratio just 
estimated appears to be ~ 50% higher than the average value 
in the diffuse ISM, X(H)/E(B - V) = 58 x 10“ cm"^ mag"' 
(Bohlin et al. 1 978). and a factor of 3 higher than what is inferred 
for the dark phase in local clouds from y-rav measurements, 
V(H)/E(B - V) == 30 X 10“ cm"2 mag"' (Grenier et al, 2005; 
Abdo et al. 2010f). The discrepancy is confirmed by the extinc- 


The same level of uncertainty would affect the Xco derivation from 
another total gas tracer .such as the dust cofumn-demsity. 
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Fig. 13. Top: emissivity per Av,esc unit versus emissivity per hy- 
drogen atom in the Cygnus complex (for Ts = 250 K). The 
points correspond to the different energy bins; the emissivities 
decrease with increasing energy. The green line gives the best 
linear fit taking into account uncertainties on both axes. Bottom: 
residuals in units of standard deviations as a function of H l emis- 
sivity. 


tion data. By fitting the latter with the H i and CO maps (§ 2.2.2) 
we obtained iV(Hl)/Av = 29.6 ± 0, 1 x 10^° cm^^ mag“* (statis- 
tical error only for Ts = 250 K). A possible explanation is pro- 
vided by an anomalous Ry ratio driven by a peculiar distribution 
of dust grain sizes (Cardelli et al. 1989), Straizys et al. (1999) in- 
deed reported an anomalous extinction law in the Cygnus region, 
showing stronger extinction in the violet and near U V region. 

The chemical state of the dark neutral gas cannot be deduced 
from y-ray observations. Whereas there are compelling theoret- 
ical and observational reasons to believe that CO-quiet H 2 is 
ubiquitous in the ISM (e.g. Wolflre et al. 2010; Magnani et al. 
2003; Langeretal. 2010), we cannot exclude that part of the 
dark gas traced by Ay excesses is missing cold atomic gas, es- 
pecially since the dark neutral phase appears at the interface be- 
tween the atomic and CO-bright phases in the nearby clouds 
(Grenier et al. 2005). Temperatures as low as 40 - 70 K were 
measured in cold H i clouds (Heiles & Troland 2003) and self 
absorption can be large when cold clouds are seen against more 
diffuse warm Hi. The Av excesses in Fig. 5 partially overlap 
an H I self-absorption feature associated with the Cygnus com- 
plex (Gibson et al. 2005, Fig. Id). Yet, the H I to H 2 transition is 
very dynamical, both in space and time, and it is difficult at this 
stage to conclude on the exact mix of cold dense H i and diffuse 
CO-quiet H 2 that forms the dark neutral phase in the outskirts of 
CO-bright molecular clouds. 

Regardless of its nature, the mass of the dark neutral gas in 
the Cygnus complex amounts at 1.4 kpc to 0.9 J x 1O*M0. 
Given an atomic mass of 5 x 10'’ Mo and including the CO- 


bright mass estimated above the total interstellar mass of the 
Cygnus complex amounts to 8);^ x 10* Mo. 

Assuming that all the dark neutral gas is molecular, we 
can calculate the molecular dark-gas fraction /do = - 

Afco)/Mnoo which amounts to'* 0.27 ± 0.02, in excellent agree- 
ment with the model by Wolfire et al. (2010). The dark-gas frac- 
tion is also consistent with that by Abdoetal. (201 Of) for the 
nearby Cepheus and Cassiopeia clouds, which have a factor of 2 
lower column densities and masses < 2% of that contained in 
the Cygnus complex. This also agrees with the prediction by 
Wolfire et al. (2010) that the dark-gas fraction is rather indepen- 
dent from the mean cloud column density and total mass for gi- 
ant molecular clouds. 


5. Conclusions 

We have performed a global analysis of y-ray emission from 
the Cygnus region measured by the Fermi LAT in the energy 
range 100 MeV-100 GeV. We built a global model for the region 
able to satisfactorily reproduce the LAT data. The model in- 
cludes extended sources that have been detected over the global 
interstellar emission model described here in association with 
the Cygnus Loop and y Cygni supernova remnants and with a 
cocoon of freshly-accelerated CRs in the innermost part of the 
Cygnus X region. They are discussed in detail in companion pa- 
pers. 

We measured the average Xco = A(H 2 )/Wco 
factor for clouds in Cygnus, finding a value [1.68 ± 
0.05 (stat.) +|5*g(Hi opacity)] x 10’° cm“^ (K km s“')“' 
well consistent with other LAT measurements for cloud 
complexes in the Local and Perseus spiral arms (Abdoetal. 
2010f; Ackermannetal. 2011). These Xco ratios, averaged 
over complexes, are, however, significantly larger than the 
values found at higher sampling resolution in nearby clouds of 
the Gould Belt (Abdo et al. 2010f). Thanks to the correlation 
between dust and y-ray emission excesses, we detected the 
presence of conspicuous masses of dark neutral gas not traced 
by the combination of the Hi and CO lines, with total mass 
~ 40% of the mass of the clouds traced by CO. The good 
correlation over three decades in energy between the y-ray 
emissivity per Ay excess unit and per H atom strengthens the 
interpretation of such excesses as due to the presence of dark 
neutral gas. The neutral gas in the Cygnus complex, combining 
atomic, CO-bright and dark masses, amounts to 8);* x 10* Mo 
at a distance of 1 .4 kpc. 

The emissivity of atomic gas measured over the whole 
Cygnus complex is consistent with that in the local interstellar 
space. We do not find evidence for the possible exclusion of CRs 
by enhanced magnetic fields in the dense clouds. The emissiv- 
ity per hydrogen atom compares with LAT estimates in other 
regions of the local and outer Galaxy, regardless of differences 
in gas surface density by about one order of magnitude and in 
Galactocentric radius by ~ 6 kpc. This uniformity does not sup- 
port models based on the dynamical coupling of CRs with mat- 
ter densities or predicting a strong emissivity gradient toward the 
outer Galaxy. 

The CR population averaged over the scale of the whole 
Cygnus complex (~ 400 pc) is similar to the Local Spur av- 
erage, in spite of the embedded regions of conspicuous massive- 
star formation and potential CR accelerators. Their impact on the 

The dark neutral gas fraction is very stable against the choice of H I 
spin temperature, therefore the error on the dark-gas fraction is statisti- 
cal only. 





CR population is detected only in the innermost region bounded 
by the ionization fronts from the massive stellar dusters over a 
scale < 100 pc (Fermi LAT collaboration, submitted). No coun- 
terpart to the broadly distributed excess of y-ray emission seen 
at energies > 10 TeV at 65° < / < 85° (Abdo et al. 2007, 2008) 
is detected at GeV energies so far. 
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